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1. Introduction 

The minimal supersymmetric standard model (MSSM) consists of a supersymmetric 
extension of the standard model, with the addition of a number of dimension 2 and di- 
mension 3 supersymmetry-breaking mass and interaction terms. It is well known that the 
MSSM is not, in fact, the most general renormalisable field theory consistent with the 
requirements of gauge invariance and naturalness; the unbroken theory is augmented by 
a discrete symmetry (.R-parity) to forbid a set of baryon- number and lepton- number vio- 
lating interactions, and the supersymmetry-breaking sector omits both .R-parity violating 
soft terms and a set of "non-standard" (NS) soft breaking terms. There is a large litera- 
ture on the effect of R-parity violation; a recent analysis (with "standard" soft-breaking 
terms) and references appears in Ref. for earlier relevant work see in particular ||. 
The need to consider NS terms in a model-independent analysis was stressed in Ref. ||; 
for a discussion of the NS terms both in general and in the MSSM context see Ref. 
however in this paper we shall ignore the NS terms. 

The unification of the three gauge couplings in the MSSM at a scale of around Mx ~ 
10 16 GeV provides compelling evidence both for super symmetry and for the existence of 
an underlying unified theory. We shall consider the standard mSUGRA scenario where 
we assume just three parameters at the unification scale, namely universal scalar and 
gaugino masses, and a universal trilinear scalar coupling, mo, an d A respectively. The 
remaining parameters are tan/3 and sgn(p). The complete mass spectrum is determined 
by running the couplings and masses from Mx to Mz (taking the quark masses and gauge 
couplings at Mz as additional inputs). There is an extensive literature on this process in the 
R-parity conserving (RPC) case and some pioneering work in the R-parity violating (RPV) 
case. In particular, Ref. |1J contained a comprehensive analysis of RPV effects on various 
scenarios, using full one-loop /3-functions for RPV parameters, and additionally including 
2-loop RPC corrections for RPC parameters. Our purpose here is to make available the 
full 2-loop /3-functions for both RPC and RPV parameters and to explore the effect of 
incorporating these full /3-functions on a representative sample of the scenarios considered 
in Ref. ffl; in particular neutrino masses and the nature of the lightest supersymmetric 
particle (LSP). We shall not actually present the /^-functions explicitly but rather refer the 
reader to a website where they can be accessed for the most general case (including a 
general 3x3 matrix of Yukawa couplings). 
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2. The Soft /3-functions 

For a general N = 1 supersymmetric gauge theory with superpotential 



W {cj>) = \^&<f>j + lYt^fafafa, 
the standard soft supersymmetry-breaking scalar terms follows 

y soft = (¥ ij + i^Mjfa + cc.) + (m 2 )%0i^', 



(2.1) 



(2-2) 



where we denote cj) 1 = etc. 

The complete exact results for the soft /3- functions are given by[lC]] — []T2 
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where 7 is the matter multiplet anomalous dimension, and 
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Here M is the gaugino mass and Y ljk = (m 2 ) l iY jkl + (m 2 ) j iY lkl + (m 2 ) k iY tjl . Eq. flj^) 
holds in a class of renormalisation schemes that includes DRED'[]T3[], which we will use 
throughout. Finally the X function above is given (in the NSVZ scheme |JL4| ) by 

^ 3 S 



X NSVZ = _ 2 9 



where 



16tt 2 [l-2g 2 C(G)(16ir 2 )- 1 ] 
S = r-hr[m 2 C(R)} - MM*C(G), 



(2.5) 



(2.6) 



C(R),C(G) being the quadratic Casimirs for the matter and adjoint representations re- 
spectively. There is no corresponding exact form for X in the DRED' scheme [fT3| ; however 
we only require here the leading contribution which is the same in both schemes; the sub- 
leading DRED' contribution is given in Ref. [ff5| . These formulae can readily be specialised 
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to the case of the RPV MSSM and their implementation can be automated; in our case 
we used the FORM package. (We have also implemented this procedure up to three loops 
for the RPC MSSM []l6| , and made the results available on another website fi% .) 

In our analysis we also include "tadpole" contributions, corresponding to renormalisa- 
tion of the Fayet-Iliopoulos (FI) D-term at one and two loops. These contributions are not 
expressible exactly in terms of (3 9i , 7; for a discussion see Ref. [|TS|. For universal boundary 
conditions, the FI term is very small at low energies if it is zero at gauge unification. 



3. The R-parity Violating MSSM 

The unbroken Af = 1 theory is defined by the superpotential 



where 



and 



W = W 1 +W 2 , 



W x = Y u Qu c H 2 + Y d Qd c H 1 + Y e Le c H x + plH x H 2 



(3.1) 

(3.2) 
(3.3) 



W 2 = \(A E )e c LL + \{Ku)u c d c d c + (A D )d c LQ + k^H^ 

In these equations, generation (i,j ■■■)■, SU2(a,b ■ ■ ■), and SUs(I,J---) indices are con 
tracted in "natural" fashion from left to right, thus for example 



A D d c LQ = e a6 (A D )^V)^ a Q 6 /- 



(3.4) 



For the generation indices we indicate complex conjugation by lowering the indices, thus 

pQ y = (y tt *)«. 

We now add soft-breaking terms as follows: 



L i = J2 m l ( t ) * ( t' + 



\H X H 2 + ^MiXiXi + h.c. 



i=i 



+ [h u Qu c H 2 + h d Qd c H x + h e Le c H x + h.c] , 
L 2 = m 2 R HlL + m 2 K LH 2 + \h E e c LL + \h v u c d c d c + h D d c LQ + h.c. 



(3.5) 



We shall also use the notation 



A ljfc = {K E ) 3 , X' ljk = (Ad) 3 , X" ljk = (Au) 



(3.6) 



with h, h! and h" denned similarly in terms of Jie, ho and hjj respectively. Note that 

^jik _ _^ijk y/ikj _ _yijk ^3 ^ 

with similar symmetry properties for h and h" . 

It can be convenient to define £q =0 ...3 = {Hf, ££=1,2,3}- The couplings A a/ 3fc, A^ Q 
are then defined so as to subsume Ag, Y e and A^, respectively; i.e A^o/c — ^eiit, Kjo = 
—Ydij- hapk, h' i j a are defined similarly. In the same spirit we define [i a = and 
b a = {ml, (nix).}; and finally m 2 c incorporates m 2 L , m 2 R and m 2 Hi . 



4. RGE Running and the Mass Spectrum 

The DR dimensionless couplings at M z are determined from the MS gauge couplings 
and the physical quark masses by incorporating supersymmetric threshold corrections. The 
boundary conditions on the soft parameters and masses are imposed at the unification scale 
M X - As mentioned earlier we adopt mSUGRA boundary conditions at M x , so we take 

m {M x ) = rriu(Mx) =mj(M x ) = m L {M x ) = m~ e {M x ) = m l, 

(4.1) 

m Hl =m H2 = m , 
where 1 is the 3x3 unit matrix in flavour space. 

Ki(M x ) ={m 2 R UM x ) = {m 2 K UM x ) = 0, 

(4.2) 

Mi(Mx) =M 2 (M X ) = M 3 (M X ) = mi . 
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Finally we define 

h u (M x ) =A Y U (M X ), h d (M x ) = A Y d (M x ), h e (M x ) = A Y e (M x ), 

(4.3) 

hu(M x ) =A Au(M x ), h D (M x ) = A A D (M X ), h E (M x ) = A A E (M X ). 

After running all the couplings from M x to Mz, the sparticle spectrum can be com- 
puted. Because of the interdependence of the boundary conditions at Mz and M x (the 
threshold corrections depend on the sparticle spectrum; the unification scale depends on 
the dimensionless couplings) we determine the couplings by an iterative process, reimpos- 
ing the respective boundary conditions at each iteration. We define gauge unification to 
be the scale where a.\ and «2 meet; we speed up the determination of this by (at each 
iteration) adjusting the unification scale using the solution of the one- loop (3- functions for 
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the gauge couplings from the previous value of the scale. We employ one-loop radiative 
corrections as detailed in Ref . . A particular subtlety in the RP V case is that the RGE 
evolution of k depends on jj, and that of m 2 K on \x and B. Therefore it is not sufficient (as 
in the RPC case) to determine \x{Mz) and B{Mz) after the iteration, from the electroweak 
breaking conditions; rather, [i and B must be included in the iteration process to establish 
values of /i(Mx) and B(Mx) which are compatible with the other boundary conditions. A 
second complication in the RPV case is the possibility of sneutrino vevs Vi, which satisfy 



v 2 = v 



1 9l 



(4.4) 



where Vd, u are the H\,2 vevs, tan/3 is defined as usual to be 



tan (3 = — 

Vd 



(4.5) 



and with our conventions v = 174GeV. Then at each iteration, \x{Mz) and B(Mz) are 
determined fromH 
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(4.7) 
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with AV being the one-loop corrections to the scalar potential (we assume the sneutrino 
vevs are real). Next the sneutrino vevs may be determined from 



(MDijVj = - [(m 2 R )i + fi*Ki] v d + (m 2 K )iV u - - 
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(4.f 



(4.9) 



Here g is the U\ electroweak coupling (usually written g'). The one-loop corrections to the 
effective potential for the RPV MSSM which appear in d ^ were obtained from Ref. 



We have included the squark contributions from Ref. [p0 [, correcting an obvious typo (a 
missing "In"); the next most significant corrections, from charged slepton/Higgs, given 
there seem clearly wrong on dimensional grounds and we have omitted them; they are 
much smaller in any case. If (as we do in the neutrino mass calculation) we impose 
electroweak symmetry breaking at the supersymmetry scale -Msusy (defined here as the 
geometric mean of the stop masses) then the effect even of the squark contributions from 
Ref. [5t| is negligible. For g^ v we have used the RPC corrections given in Ref. ||19|| . (For 
the calculations of selectron, stau and stop masses given later we incorporate one-loop 
threshold corrections and therefore the choice of EWSB scale should be less significant; 
and in fact we choose to evaluate the sparticle masses at their own scale.) 

Our philosophy throughout is to investigate qualitative effects, particularly of using 
two-loop rather than one-loop /3-functions. Therefore we have made various simplifications 
in our procedures. ADD consider three standard forms for the relation between the weak- 
current and quark-mass bases for the couplings, where there is either no mixing, or the 
mixing is all in the down-quark sector, or all in the up-quark sector. We have assumed 
that the Yukawa matrices are diagonal in the weak-current basis both at the GUT scale 
and at the weak scale. This corresponds to assuming a trivial CKM matrix, Vqkm = 1 
at the weak scale. We are also neglecting the generation of off-diagonal Yukawa couplings 
in the evolution from Mz to Mx (an effect which we believe is negligible to the accuracy 
at which we are working). 



5. Neutrino Masses 



Here we set bounds on the couplings A, A' from the cosmological neutrino bound. 



Combining the 2dFGRS data|2l|] with the WMAP measurement [p2| one gets a bound on 
the neutrino mass 

^2m Ui < 0.71eV. (5.1) 



The neutrino mass is given by 



m,. 



fj,(M ig Z + M 2 g 2 ) E - = i A 



2 {v u v d {M l9 l + M 2 g 2 ) - vM 1 M 2 ) 



(5.2) 



where 

Ai = Vi-v d — ■ (5.3) 
I 1 

A single non-zero RPV coupling at Mx will generate non-zero k, m? R and m 2 K leading to a 
non-zero neutrino mass. We follow Ref. |1]] in choosing the SPSla mSUGRA point, which 
has the following parameter values at Mx'- 

m = lOOGeV, mi = 250GeV, A = -lOOGeV tan/3 = 10, sign(//) = +. (5.4) 
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Eq. (|5.1|) then leads to an upper bound on the given RPV coupling. We assume that only 
one out of the set of couplings 

S\ = {A333, A 322 , A 311 , A233, A232, A131} (5.5) 

is non-zero at M z , and that only these couplings are non-zero in the running; these very 
nearly form a closed set in any case, since the only additional couplings which could be 
generated (at one loop) are A 211 , A 222 an d A121. Looking at the form of the /3-functions 
one can see that these couplings could not in any case be generated at a level close to 
their limiting values, since the coupling in S\ responsible for generating them has a much 
smaller limiting value and is additionally suppressed by small (1st or 2nd generation) RPC 
Yukawa couplings. Moreover (if we start with just one of them non-zero) these couplings 
do not generate any off-diagonal contributions to Y Uj d, e so our assumption about the form 
of these matrices at Mx is justified. 

The bounds on these couplings are shown in Table 1. 



Coupling 


1 loop 


2 loop RPC 


2 loop 




1.0 x 10" 5 


8.7 x 10" 6 


8.4 x 10" 6 


A 322 (M Z ) 
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3.2 x 10" 4 
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A 233 (M Z ) 
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5.3 x 10" 5 


5.4 x 10" 5 


A 232 (M Z ) 


1.1 x 10" 3 


1.0 x 10" 3 


9.2 x 10" 4 


Ai3i(M z ) 


2.2 x 10" 1 


1.9 x 10" 1 


1.8 x 10" 1 



Table 1: Upper bounds on A(M Z ), X'(M Z ) 
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The "2 loop RPC" column corresponds to the procedure followed in Ref. [II] , where the 
full one-loop /3-functions were used and also the two-loop RPC corrections were included 
in the /^-functions for the RPC couplings and masses. Our results (in the 2-loop RPC case) 
agree well with those of ADD, particularly for the A' limits where we agree to better than 
2%. 



6. The Nature of the LSP 

In the RPV case the LSP is no longer stable and therefore no longer subject to 
cosmological constraints on stable relics. Also the LSP need not be electrically and colour 
neutral. Once again we follow Ref. in taking for this analysis the case of "no-scale" 
supergravity, which corresponds to taking Aq = mo = 0. We shall consider the variation 
of the nature of the LSP with A231. In this case the LSP is either a stau or a selectron. 
The computation of selectron masses is in general more complex than the RPC case, since 
the charged Higgs mix with the charged sleptons, giving mass terms of the form 

fht\ 

C ch = -( }q e Ll e Rk ) M c 2 h I <g J , (6.1) 

where M 2 ^ is an 8 x 8 matrix given by 

/(m 2 )u + D b* s + D s \palH* v P 

M ch = b l + D j (™> 2 )s-f + XcrylXpSlVaVp + D l5 h ai iv a - \ ai in* a v u 



^pakVatVp Ksk V » ~ * Kjlt + K(3k^ al lVpV 7 + D 



Ik 



(6.2) 



where 



I 2\ — 2 1 1 i2 



2 

a ' 



(6.3) 



D5 =\glv u vs, 

a 

However if A231 is the only RPV coupling, the matrix is still diagonal except for the 
standard stau mixing. We have included the one-loop corrections to the slepton masses as 
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given in Ref. . Of course these omit any corrections from RP V couplings but presumably 
these will be extremely small. 

In Fig. 1 we show the variation of the nature of the LSP with tan/3 and A23i(Mx). 
Here we have used the two-loop RG evolution equations but in fact the results using one- 
loop evolution are almost identical. Moreover it is easy to check that (at least at one loop) 
if A231 is the only non-zero RPV coupling at Mz then it will remain so at all scales, so we 
can use a simplified set of /^-functions in which we only retain A231. 



0.45 



Boundary between selectron and stau LSP 



0.35 



0.25 




0.15 



0.05 



Fig. 1: The variation of the nature of the LSP (stau LSP below the line, 
selectron LSP above). 

Once our results agree pretty with those of Ref. JJ] , although our demarcation line is 
slightly lower, particularly for larger values of tan (3. 



7. The stop masses 

The bounds on the A" couplings are much weaker than for the A and A' couplings 
and in general are only set by perturbativity of the top Yukawa coupling. The situation 
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changes if a particular form is assumed for the quark mixing, such as mixing only in the 
up-quark or only in the down-quark sector. The bounds are particularly stringent in the 
down-quark mixing case. Although, as described earlier, we assume the no-mixing case, 
we expect our results to be qualitatively valid in the general case and therefore we shall 
display our results up to the perturbativity bound. We consider the dependence of the stop 
masses on A323. (In the no- mixing case it is clearly consistent to consider a single non-zero 
coupling at all scales.) The mass matrix for up- type quarks has no explicit dependence on 
the RPV couplings and so the dependence on A323 is purely an implicit effect due to the 
RG evolution. The stop masses are very sensitive to the value of the top mass; here as 
elsewhere in the paper we take m top = 174.3GeV. We see that the variation of the stop 
masses, especially the light one, on A323 is considerable. 



Stop masses versus X, 
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Fig. 2: Variation of the light stop mass with A323. 
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8. Conclusions 

We have analysed the effect of including the full set of two-loop /5-functions for Im- 
parity violating couplings in a variety of scenarios. Typically we find little difference 
between the effect of using the full /3-functions and that of using the one-loop /3-functions 
plus two-loop RPC corrections for RPC parameters; though as we see in Table 1 there is 
quite a substantial difference between the bounds on RPV couplings obtained using the full 
two-loop /^-functions and those obtained using the full one-loop /5-functions-and of course 
it is desirable from the point of view of consistency to use the full set of /3-functions. In any 
event, we hope that future analysts will find the availability of the full set of j3- functions 
for the most general R-parity violating version of the MSSM to be a useful resource^]. 
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